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THERMALLY INDUCED OSCILLATIONS IN 
CRYOGENIC SYSTEMS* 


R. S. Collier 


This work was done under NASA Order W 13, 300 for the period 
November 1, 1970 to October 31, 1971 as a part of Task V- Slush Hydro- 
gen Instrumentation. 

One device — a carbon film sensor — has been found to possess a 
peculiar feature in addition to its demonstrated capability to discriminate 
among vapor, liquid, and slush phases. When a carbon film resistor is 
placed in an environment of slush hydrogen in a dc power mode, fluctu- 
ations of the carbon film temperature are detected which are unique to 
the slush environment. 

These temperature oscillations can be correlated with fluctuations 
previously observed in research on the characterization of slush. This 
characteristic of slush hydrogen - the susceptability of the mixture under 
certain conditions to exhibit violent pressure excursions in valve stems, 
pressure gage lines,- etc. — is probably the most singular difference be- 
tween the technologies of slush and liquid hydrogen. 


The similarity of the response of the carbon film sensor to ther- 
mally induced oscillations previously noted was suggested, and a study, 

... both experimental and analytical;' was proposed to: 

1) Determine the nature of the temperature oscillations ob- 
served on the carbon thin films; 

2) Proceed to perfect this point sensor as a slush hydrogen 
detector if these oscillations can be shown to have a firm physical and 
theoretical basis; 

♦Work carried out at the National Bureau of Standardr under the sponsorship 
of National Aeronautics and Space Administration (SNSO) Order No. W13,300. 


3) Determine if there it a connection between the two observed 
phenomena — the temperature fluctuations of the carbon film and the so- 
called "thermal oscillations" observed in slush studies; 

4) Assess whether these phenomena are potential problem areas 
in the extension of the use of slush hydrogen as a space fuel and propellant. 

The first three goals have been pursued during this contracting period 
and the final results of this effort are contained in three separate papers 
contained in this report; the papers are being individually submitted to 
technical journals for publication and are titled as follows: 

Part I Temperature Oscillations of a Heated Surface in Slush 
Hydrogen. 

Part II The Expansion Equations of Non- Equilibrium Thermo- 
dynamics. 

Part III Surface Damping of Thermally Induced Oscillations in 
-Cryogenic Systems. ....... 
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Part I 


Temperature Oscillations of a Heated Surface 
in Slush Hydrogen 

R. S. Collier 
Cryogenics Division 

Institute for Basic Standards, National Bureau of Standards 
Boulder, Colorado 80302 

ABSTRACT 

\ 

_ It was found that the surface temperature of a steadily heated 

cylinder oscillates with an amplitude of about 0. 25 K and an oscillation 
frequency of about 8 Hz when immersed in slush hydrogen. The ampli- 
tude and the period of oscillation are relatively independent of heat flux 
at the surface in the range 0, 25 - 1. 5 watts /cm 3 . The experimental 
results are consistent with a theoretical model which predicts oscillations 
in the thermal expansion of a liquid film located between the heated sur- - 
face and the slush. . 

'n 

Key Words: . Carbon thin film; slush hydrogen; thermal oscillations. . 
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Introduction 

A pyrex cylinder, approximately 1 cm long and 2 mm diameter, 
was coated with a vacuum deposited carbon film about 500 A thick. The 
resistance of the carbon film was approximately 2300 Cl in slush hydrogen 
and was measured by attaching electrical leads to each end of the cylinder. 
The temperature sensitivity of the film in the liquid hydrogen range was 
approximately 50 fl per kelvin. The carbon film was used as a constant 
current heat source with surface heat flux between 0. 25 watt/cm 8 and 
1. 5 watt/cm a . When this heat source was immersed in slush hydrogen, 
there was an 8 - 10 Hz voltage oscillation corresponding to a temperature 
oscillation about 0. 25 K in amplitude. This temperature oscillation is 
consistent with a theoretical model of a "relatively incompressible, slowly 
accelerating fluid" which predicts oscillations in the thermal expansion 
of a fluid near a steadily heated surface^ When the heat source was 
immersed in liquid there were occasional oscillations at about 16 Hz and 

the amplitude of these oscillations was much smaller. 

/ 

Experimental ' 

The carbon film used in this experiment is similar to those 
developed for fluid phase discrimination^’^, where the film is subjected 
to a constant current pulse and the transient response is analyzed to de- 
termine whether or not the sensor is in the liquid or gas phase. When 
the sensor was immersed in slush hydrogen, the transient response as 
it appears on an oscilloscope trace is shown in figures 1-3 for three dif- 
ferent power levels. The initial part of the transient is the same in the 
liquid and in the slush; indicating that the sensor is surrounded by a layer 
of liquid; the surface temperature then becomes slightly warmer in the 
slush than in the liquid and steady state oscillations appear which persist 
indefinitely. 
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Theoretical Interpretation 

It has been shown that these steady state oscillations can be 
described as approximate solutions to the conservation equations of non- 
equilibrium thermodynamics ^ 1. The period of oscillations, t, is given 


10. 2d 3 

T * (1> 
m 

where a is the thermal diffusivity, P is the Bulk Coefficient of Thermal 
Expansion, AT is the peak to peak surface temperature oscillation, and 
d is the thickness of the liquid film on the surface of the sensor. The 
thickness, d, is the only quantity which is not directly measurable. 

For liquid hydrogen at 14 K, op = 2 x 10' 4 cm a /sec K. It is seen 

in figures 1-3 that AT (rs 0. 2 - 0. 3 K) is relatively independent of the 

m 

power input which theoretically predicts, by equation (1), that the period" - 
of oscillation, t, is relatively independent of power input. This is seen 
to be true, 1 /t es 8 - 10 Hz for all three power levels. From these numbers 
the liquid film thickness is calculated to be d w 0. 8 X 10" 3 cm. 
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Part II 

The Expansion Equations of Non- Equilibrium Thermodynamics 

R. S. Collier 
Cryogenics Division 

Institute for Basic Standards, National Bureau of Standards 
Boulder, Colorado 80302 





ABSTRACT 

Definitions of "Relatively Incompressible" and "Slowly Accelerating" 
are used to convert the differential forms of the conservation equations to a 
pair of coupled non-linear partial differential equations in the density, p, 
and the thermal expansion velocity, v. An approximate solution to these 
"expansion equations" indicates that oscillatory solutions can exist for the 
thermal expansion of a fluid near a steadily heated surface. This example 
suggests several conjectures which are posed subject to a more exact 
treatment of the expansion equations. 

Key words: Expansion equations; non-equilibrium thermodynamics; 
slush mixtures; thermal oscillations. 
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Introduction 


Six independent fluid variables of non-equilibr’um thermodynamics 
are given by the density, p, the temperature, T, the pressure, P, and 
the velocity, v, which is a three dimensional vector. These variables 
are related by the conservation equations of mass, momentum, and energy 
and by the equation of state. The conservation equations have several 
forms which are equivalent under certain conditions. The differential 
forms, neglecting viscous and body forces, are given in rectangular 
coordinates as follows^: 


Conservation of Mass: 


5p 

at + 7 * pv = 0 


or 


Dp + p7»v = 0. 

B. Conservation of Momentum: 

pDv = -VP 

or 


5t ( Pv) + 7 (tt) = 0. 


C. Conservation of Energy: 


pC DT = pQ + 7«k7T - T(~) 7 «v 


or 


.ir^) 


pC p DT = pQ + 7«k7T - - DP. 


( 1 ) 

( 2 ) 

(3) 

(4) 


(5) 


( 6 ) 


In these equations, 7 is the vector differential operator, D is the 
convective differential operator given by D = 3/St + v»7, C and C 


are 


the specific heats of the fluid at constant volume (density) and constant 

pressure, respectively. Q represents the heat input per unit mass, and 

n is the momentum flux density tensor which is given in this case by 

(n).. = P6. . + pv.v,; here, v. and v. are velocity components and 6.. is 
ij ij i j 1 J r ij 

the Kronecher Delta. The viscous and body forces are neglected for 
simplicity and do not have a material effect on the arguments presented 
here, except that it will normally be much harder to obtain solutions to 
the resulting equations if these forces are included.. - 


A Relatively Incompressible Fluid 
The equation of state for small density excursions about p^ can 


be written as 


where 


p = p Q - pB(T - T o > + pK(P - P q ) 


8 ■ p^dT 


H is the Bulk Coefficient of Thermal Expansion 

P ' 


K = i* the Isothermal Compressibility Coefficient. 

P T 


In the remainder of this discussion, the situation is considered 
where the lasttermin equation (7) can be neglected; that is, to a good 
approximation, p and its derivatives can be given by 

» = p„ - »8(t . T d ) ) ; ' 


B • - pelf, pp = - oBPT 


where p , and T . are constants. In this case, the fluid is considered 
d d 

incompressible relative to its thermal expansion. This situation is of 
considerable importance in many applications involving gases and cryo> 
genic fluids where the heat transfer is high, $_ is large, and the pressure 
effects are small. 1 ...... . 






Under equations {8), the energy equation can be easily simplified. 
Combining equations (5) and (2), it follows that 


pC p DT ‘ ^ ='pQ + 7.krr, ( 9 ) 

' p p 

and using (8), equation (9) becomes 

(p C o + T (f)0 DT = pQ + 7%k7T - (io) 

p p 

The coefficient of DT in equation (10) is equal to pC_. This follows from 
the basic thermodynamic identities^ 



Hence, the energy equation becomes 


pC p DT = pQ + 7«k7T. (13) 

It is also possible to eliminate DT and 7T from this equation by using 
(8) and (2), and it follows that 



7»v + 7»^7p. 


(14) 


This appears to be the most useful form of the energy equation. 
The physical content of equation (14) is that the power input density at a 
point is divided between the thermal conduction output and the convective 
output due to thermal expansion. 


It is seen by comparing equations (6) and (13) that the hypothesis 
of "relative incompressibility, " equations (8), implies DP = 0 or 


dP 

at 


— + v. VP = o. 


(15) 


The final equation to complete this set is obtained by operating on 
equation (1) by d/at, switching the time and space derivatives in the second 
term; also operating on equation (4) by V*. These two results may be 
combined to give 

a s p 

- V. V(TT ) = 0 (16) 

where 


v. v(tt) = V 2 P+.E 


^ dx 1 dx^ 


pv.v.. 
i J 


Equation (16) is normally used to predict non-linear effects associ- 
ated with sound waves; it is relevant to note that equation ()6) predicts 
the existence of shock waves, but does not describe themf 3 !. 

The relative incompressibility hypothesis thus implies the follow- 
ing equations for non-equilibrium thermodynamics: 

PC P k 

pQ - — — V»v + V- — VP 
. . . B dS 


- V.V (tt) = 0 


(17) 


DP 0 or. — - + v-VP = 0. 

o t 


A Slowly Accelerating Fluid 


The pressure P can be eliminated from these equations if we are 
willing to concentrate on the dynamics of the thermal expansion process 
near a heat transfer surface. In this situation the expansion is controlled 
by thermal diffusion which is a relatively slow process as compared to the 
propagation of sound waves. Hence, we can neglect the effect of sound 
waves by assuming that the fluid is slowly accelerating and that the pres- 
sure and velocity are related by Bernoulli's equation 

P = P - ipv a (18) 


fcaS.ilkv' 


where P^ is a constant ullage pressure. 

In this case, the first two non- equilibrium equations (17) become 
a pair of coupled non-linear equations in p and v 

’ '.V-- -pc: 


P Q = 7.v+7.^7p 


at 


3 - 7*7(TT) = o 


(19) 


where 


7» 7(tt ) = - ^7 a (pv 3 ) + 2 


ij Sx ax 


pv.v., 

i j 


•_ Equations (19) will be referred to as "the expansion equations", and 

it will be shown by an example that the thermal expansion described by 
these equations can have oscillatory, wave-like properties near a heat 
transfer surface. 

Thermal Oscillations - 

Although sound waves have been neglected, they have not been 
forgotten. Indeed, an important effect is the coupling of the thermal 


I 





expansion waves to the possible acoustic modes of the system. Thes'' 
are the effects known as thermo-acoustic oscillations or thermally induced 
pressure oscillations^. In describing this coupling, both hypotheses 
made here will have to be relaxed far enough to include, at least, linear 
sound propagation. This discussion is relevant to a consideration of these 
effects, however, since the thermal expansion considered here is a possi- 
ble driving mechanism for the thermo-acoustic oscillations. 

Example - Heat Transfer in Slush 


As an example of an approximate integration of the expansion 
equations, consider a slab (film) of liquid between a flat heater at x = 0 
with heat flux q(t) and a porous surface at x = d kept at a constant tem- 
perature T . . In this case the expansion equations for 0 < x < d become 


_£ -Sv + JL = o 

3 3x dx p3 3x 



( 20 ) 

(2D 


where v is described by a single component in the x direction. The 
boundary conditions are given by 


v 

p(d, t) 


The quantities 
equation (2) becomes 



7b f£ ‘ <l " 1 *•* = 0 ■' < 22 > 

(a constant) at x = d. (23) 


K - - • 

and — — will be assumed constant so that 
PP 


where O = —p — is the Thermal Diffusivity. Equation (21) has a first 
pC P 

integral given by 

v + ||£= f (t, (25) 

where f is an arbitrary function of t. The application of the boundary 
conditions (22) to equation (25) gives 


v _ Pqft) . 

pc p p ax • 


(26) 


If the spacial derivatives of v 3 are neglected, the system of equations 

become 


v _ Bdft) er Sp 

pc p ' p ax 


d a P . 

at 3 


-iv a 



(27) 



This is a quasi-linear, coupled, set of equations and there exist separable 
solutions of these equations which satisfy the boundary conditions. We 
look for solutions of the form 

;o(x,t) = ■?» Q (t) + - R(t) S(x) (29) 

p p ° P c p 


where p Q (tj is the density at x = 0, which is deterTnined by measuring the 

surface temperature and using the relation (8); R and S are functions to 

be determined. It is convenient to consider p (t) as the sum of two 

o 

quantities 

P o (t) = P oo (i) + Ap o (t) (3°) 

where p QQ (t) depends only on the nature of the heat source, q(t), and 


Ap Q (t) is the density (and hence temperature) change caused by the dynamics 
of the expansion process. If p oQ (t) * s given by 


o a 

~p (t) a -p, _ 

p oo p d pC 


then it follows from (23) and (29) that 


Eq(t)d 

p 


(3 1) 


- Ap (t) = R(t) S(d). 
P o 


Also by substituting (29) into (27) it follows that 
v(x,t) = R(t)S'(x) 


(32) 


(33) 


where the prime, ' , indicates differentiation with respect to the argument. 


If the heat input is constant, q(t) = q , then n is just the static 
r o K oo 

equilibrium surface density. Thin case will now be considered: 

It follows that equation (28) becomes - _ ... 

R'(S(d) - S(x)) - *R 3 (S(d) - S(x))' 8 'S(d) - S(x))' = 0 (34) 

where we note in passing that (S(d) - S(x) ,z = (S'ts)) 2 . Defining AS(x) = 
S(d) - S(x), it follows that the boundary conditions for AS(x) are 


AS(0) = S(dJ = S d . AS(d) = 0. (35) 


Here we note that S (0) = 0. 


Separating variables in (34), ordinary differential equations for 
R and AS are obtained 

R* - bR s . - . 0 -- ...... - .. (36) 

(AS) ,a (AS)* - 2b (AS) = 0 (37) 

where b is a separation constant. These equations are directly integrable 
in terms of Elliptic Functions 


( 38 ) 


The first integrals of (36) and (37) are, respectively, 

R /a = fc b (R 4 + E) 

(AS) '* = 4b (( AS) s f F) (39) 

where E and F are arbitrary constants. Equations (ZZ), (33), and (35) 

imply from the conditions at x = 0 that F = - S , a , and equation (39) be- 

a 

comes 

S' 4 = 4b S(S - 2S ,). (40) 

a 

A graph of S(S - 2S.) shows that S(S - 2S.) < 0 for S between zero 
d d 

and no matter whether sign of S is positive or negative; hence, the 
reality conditions of equation (40) dictate that 0 a b H - K®, K 2 0. If b = 0, 
equations (33) and (40) indicate static equilibrium; if b < 0 then an analysis 
of equation (36) shows that the solution is oscillatory as far as the time 
development goes. In this case, equations (38) and (40) become 

'• • R ,a = i* a (R 4 - R 4 ) • (41) 

m _ • . 

S /4 = 4K a S(2S d - S) (42) 

where l v m ( x )| = I R m S *( x ) I is the maximum speed of expansion at the 
point x and serves to define the constant, R . The solution tc equation- 
(41) is given by - - - - ... 

R = R cn (R Kt, 1//Z) (43) 

mm 

where cn (R Kt, \//l) is the elliptic cosine of (R Kt) of modulus \ I fi. 
m r m 

This is similar in form to the cosine function except in this case, the 
oscillations are more sawtooth in shape. A graph of this function is given 
in figure 4. The period of oscillation, t, is given by 



4 Art J_N _ 4 x 1. 85 

R K * KZ’/lJ = R K 


where F(tt/2, 1//Z) is the Elliptic Integral of the First Kind of modulus 
1//2 and argument tt/2 (Complete Elliptic Integral). As is normal for 
non-linear oscillations, the period depends on the amplitude, R , and in 
this case is inversely proportional to R^ 

The separation constant K can be evaluated by solving (42) and 
applying the boundary conditions at x = d. In this case, it is interesting 
to do this in order to see how the period of oscillation, t, depends on the 
other parameters of the problem. 

The solution to equation (42) is given by 

±/T(Cx = /Z3T [(V(tt/ 2, 1//2) - F(cp, l//Z)J 

+ 2(E(cp, 1 IfZ) - E(n/2, 

where F and E are Elliptic Integrals of the First and Second Kind respec- 
tively, of modulus 1//2; the argument ep is given by 

rp(x) = cos 

and it is seen that cp(d) = 0. Therefore, the separation constant, K, is 
given by 

c 

tS 


tc = "jt [f(tt/2, - 1 // 2 ) - 2E(n / 2 ; \//i ) j 


( 48 ) 


From equations (44) and (47), the period, r, is given by 

„ _ 10. 2d 3 

R S, ’ 

m d 

or tn terms of the maximum dynamic density excursion at the surface 
AP om * see e< l uation (32)), 


T 


10. 2 d 3 


a 

~ Ap 
p om 


(49) 


and finally, in terms of the maximum dynamic surface temperature ex 
cursion, = 0AP /P (see equation (8)), 


T _ 10- 2 d a 

08 AT ' (50) 

m 

This expression has been experimentally verified as being approxi- 
mately correct for slush hydrogen^ 6 ]; it is seen that the oscillation period 
is also relatively independent of the amount of heat input, and this has 
also been verified f 6 1 (In a more accurate treatment of ^expansion 
equations if is likely that AT" m will have a weak dependence on ) 

Discussion 

The example outlined above, specialized as it may be, allows the 
formulation of several conjectures which should be checked by a more 
exact treatment of the expansion equations and also could form the basis 
of an experimental program: 

\ ' ; , " . ’ ' • . i 

i) Oscillations in the thermal expansion can occur in any fluid 
within a distance of the order of a thermal diffusion length of the heat 

transfer surface, provided the thermal expansion coefficient is high enough. 
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2) The frequency of the oscillation is relatively independent of 
the magnitude of the heat flux at the surface. 

3) T* is seen that the initial conditioas of the example are for- 
mulated assuming that the fluid is in a condition of steady state oscillation; 
it is quite likely that the existence of oscillations depends on the nature 

of the initial conditions; and in particular, what is q(t)? From this pc:nt 
of view it is also quite likely that mechanical perturbations can initiate 
oscillations; this has been seen experimentally^]. 

4) It is possible that the one dimensional nature of the example 
is essential for the existence of oscillations. If the expansion velocity 
has more than one component, then the tensorial properties of rr make it 
quite unlikely that separable solutions could be obtained. Although separ- 
ability is not, a priori, necessary for the existence of oscillations, it 

certainly enhances the possibility. This conjecture has been checked 
r 8 i 

experimentally 1 J by cutting machine threads into the inside of siphon 
tubes^}, and also by threading a valve stemt^®}. In these cases the os- 
- cillations are damped; possibly because of a three component expansion - 
velocity. 

5) It is likely that viscous forces, if they are large, could 
dampen the oscillations both from the point of view of separability and 
also as an energy absorbing mechanism. 

6) An analysis of equation (45) near x = 0 shows that there is 
a weak singularity in the density at the heated surface suggesting that the 
oscillations in an extended fluid may be associated either with shock waves 

— or weak shock waves. If this is the case, it may be necessary to use the 
integral form of the conservation equations in order to obtain a complete 
solution to the problem. 

- - -7) It is likely that these expansion oscillations may excite possible 

acoustic modes of a system (thermo-acoustic oscillations). In this case, ' 
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the oscillation frequency will probably be determined by the acoustic mode. 
The expansion oscillations may be at a different frequency, but since this 
oscillation is non-linear it would have a broad excitation spectrum as a 
forcing mechanism for the acoustic mode. 
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ABSTRACT 

Theoretical and experimental investigations have shown that thermal 
oscillations originating at a steadily heated surface depend on the nature 
of the thermal expansion processes at the surface as described by the 
expansion equations of non- equilibrium thermodynamics. These equations 
indicate that thermal oscillations are stronger and easier to initiate when 
the heat transfer surface has a smooth, regular geometry; and are weak 
and hard to initiate and maintain if the heat transfer surface is rough or 
irregular. Simple experiments are described which substantiate these 
findings. 
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Introduction 


The existence of thermally induced oscillations is a recurring 
problem in liquid helium and slush hydrogen cryogenic systems. These 
low frequency oscillations are normally observed as pressure fluctuations 
and high heat transfer in gauge lines, valve stem clearances, relief valves, 
siphon tubes, etc. It is possible that large scale oscillations which are 
occasionally observed in cryogenic systems may have the same origin. 

They have also been observed as surface temperature fluctuations, using 
a resistive carbon film sensor immersed in slush and liquid hydrogenl 1 !. 

These oscillations are always associated with a region of the system 
where there is heat conduction from a solid constraint into the system; 
hence, the term "thermally induced oscillations. " Various methods have 
been proposed for controlling or damping the oscillations which include 
drilling small vent holes in closed tubes teeing- in Helmholtz resonators 
to external linesfa], and inserting woolen threads into siphon tubesM. In 
this note, a method is proposed which could be applied to many different 
situations. - - - 

The conservation equations of non- equilibrium thermodynamics 
suggest that thermally induced oscillations originate as non-linear oscil- 
lations in the thermal expansion of the fluid near the heat transfer surface^. 
These non-linear oscillations can occur if 1) there is a weak singularity in 

the expansion velocity at the heat transfer surface and 2) there is a- long 

range coherence of the expansion velocity within a distance, d, from the 
surface. The distance,, d, depends on the heat transport properties of the 
fluid and the geometry of the surface; for example, for a flat surface in 
liquid hydrogen, da» 10“*cm. 

It is tempting to suggest that a surface roughness on the order of 
d will excite random components of the expansion velocity destroying the 



long range coherence and thus damping the oscillations at the source. 
Accordingly, it was proposed^! that a screw thread cut into the inside of 
a siphon tube might be of sufficient roughness to dampen the oscillations. 

Four brass tubes, 1 /4-inch OD and 5 inches long, were fabricated 
as follows: three of them were drilled with a #46 (0. 0810 inch) drill and 
two of these were tapped with a 3-56 tap (75% full thread); the fourth tube 
was drilled smooth with a #39 (0. 0995 inch) drill corresponding to the OD 
of the tap. These brass tubes were soldered to four stainless steel siphon 
tubes 1 /8-inch OD, 0. 010 inch wall, and three feet long. They were in- 
serted (after precooling) into a 25 liter liquid helium storage dewar which 
was about half full of liquid. Oscillations were detected as pressure os- 
cillations in a small cup which was soldered to the other end of the stain- 
less steel tube. For the two smooth drilled brass tubes, the oscillations 
were very strong, easy to excite, and persisted indefinitely. For the two 
tapped brass tubes, the oscillations were very weak, hard to excite, and 
occasionally would damp out completely. 

- The principles outlined above were applied to an existing problem 

involving a valve stem clearance in a slush hydrogen system [?]. The 
valve stem was a 3/32 -inch stainless rod which had a small flat Hied along 
its 4-foot length. The stem traveled through a 1/8-inch, 0. 010-inch wall 
stainless steel tube which was inserted into a 1 /4-inch stainless steel tube. 
The 1/4-inch tube was a vacuum jacket and the space between the 1/4-inch 
tube and the 1/8-inch tube was filled with helium gas at — 5 psi. The os- 
cillations were observed as vapor bubbles streaming from the L /8-inch 
tube whenever the end of. the tube was inserted into the slush, thus indi- 
cating an abnormally high heat transfer. Various efforts to open, close, 
change the gas or pressures at the other end of the tube failed to suppress 
the bubbles although the bubbling would occasionally stop after a period of 
about 20 minutes. 
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The valve stem was then threaded along its entire length with a 
3-56 die. It was found that the bubbling damped completely in less than 
one minute after inserting the 1 /8-inch tube into the slush. 
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